SPARTA (7 Jul 2014)
# flow around a sphere on a uniform grid
# particles are created initially and input at boundaries at stream velocity

variable            x index 10
variable            y index 10
variable            z index 10

variable            n equal 10*$x*$y*$z
variable            n equal 10*40*$y*$z
variable            n equal 10*40*50*$z
variable            n equal 10*40*50*50
variable            fnum equal 7.33e+21/v_n

seed	    	    12345
dimension   	    3
global              nrho 7.03E18
global              fnum ${fnum}
global              fnum 7.33e+15
timestep 	    1.E-5
global              gridcut 0.1
global              surfmax 300

boundary	    o o o
create_box          -5.0 5.0 -5.0 5.0 -5.0 5.0
Created orthogonal box = (-5 -5 -5) to (5 5 5)
create_grid         $x $y $z
create_grid         40 $y $z
create_grid         40 50 $z
create_grid         40 50 50
Created 100000 child grid cells
  parent cells = 1
  CPU time = 0.00488615 secs
  create/ghost percent = 88.0648 11.9352

balance_grid        rcb cell
Balance grid migrated 87500 cells
  CPU time = 0.038362 secs
  reassign/sort/migrate/ghost percent = 11.0297 0.88066 25.0668 63.0228

species		    ar.species Ar

mixture             air Ar frac 1.0
mixture 	    air group species
mixture             air Ar vstream 2500.0 0.0 0.0 temp 300.00

fix		    in inflow air all

collide		    vss air ar.vss

read_surf	    1 data.sphere
  602 points
  1200 triangles
  -1 1 xlo xhi
  -1 1 ylo yhi
  -1 1 zlo zhi
  0.100631 min triangle edge length
  0.00439601 min triangle area
  400 = cells with surfs
  4040 = total surfs in all grid cells
  22 = max surfs in one grid cell
  0.402525 = min surf-size/cell-size ratio
  99360 240 400 = cells outside/inside/overlapping surfs
  400 = surf cells with 1,2,etc splits
  995.862 995.862 = cell-wise and global flow volume
  CPU time = 0.022588 secs
  read/surf2grid/error/ghost/inout percent = 6.96847 9.59352 11.7805 61.6924 9.96506
surf_collide	    1 diffuse 300.0 1.0
surf_modify         collide 1 1

create_particles    air n 0
Created 955103 particles
  CPU time = 0.0248909 secs

fix                 3 balance 100 1.1 rcb part
collide_modify      vremax 100 yes

#compute             1 grid all ndensity mass u v w
#fix                 1 ave/grid 10 100 1000 c_1 ave running
#dump                1 grid 1000 tmp_flow.* id xc yc zc f_1

#compute             2 surf all press ke
#fix		    2 ave/surf 10 100 1000 c_2 ave running
#dump                2 surf 1000 tmp_surf.* id f_2

compute             g grid all n
compute             max reduce max c_g[1]

stats_style	    step cpu np nattempt ncoll nscoll c_max
stats		    100

# equilibrate flow for 1000 steps
# then benchmark for 1000 steps

run 		    1000
Memory usage per proc in Mbytes:
  particles (ave,min,max) = 12.3596 12.3596 12.3596
  grid      (ave,min,max) = 3.73848 3.73848 3.73848
  surf      (ave,min,max) = 0.0601273 0.0601273 0.0601273
  total     (ave,min,max) = 16.2536 16.2536 16.2536
Step CPU Np Natt Ncoll Nscoll max 
       0            0   955103        0        0        0           11 
     100    1.3776009   954010      287      265      112           68 
     200    2.8794088   955264      389      347      110           62 
     300    4.3786399   961823      455      403      113           72 
     400    5.9876139   976270      460      414      112           70 
     500     7.484967   991671      465      410      106           72 
     600    8.9895139   996260      498      450      119           75 
     700    10.504374   998960      490      438       97           77 
     800    12.034761   999898      490      442      123           74 
     900    13.559825  1000552      497      441      122           71 
    1000    15.076036   999719      478      428      119           69 
Loop time of 15.0761 on 8 procs for 1000 steps with 999719 particles

Particle moves    = 983950998 (984M)
Cells touched     = 1099163904 (1.1B)
Particle comms    = 2844278 (2.84M)
Boundary collides = 0 (0K)
Boundary exits    = 2735150 (2.74M)
SurfColl checks   = 63947372 (63.9M)
SurfColl occurs   = 113273 (0.113M)
Collide attempts  = 14431909 (14.4M)
Collide occurs    = 11064371 (11.1M)
Particles stuck   = 0

Particle-moves/CPUsec/proc: 8.15821e+06
Particle-moves/step: 983951
Cell-touches/particle/step: 1.11709
Particle comm iterations/step: 1
Particle fraction communicated: 0.00289067
Particle fraction colliding with boundary: 0
Particle fraction exiting boundary: 0.00277976
Surface-checks/particle/step: 0.0649904
Surface-collisions/particle/step: 0.000115121
Collision-attempts/particle/step: 0.0146673
Collisions/particle/step: 0.0112448

Move  time (%) = 7.89803 (52.3879)
Coll  time (%) = 1.99868 (13.2573)
Sort  time (%) = 4.64911 (30.8377)
Comm  time (%) = 0.199987 (1.32652)
Outpt time (%) = 0.0412275 (0.273463)
Other time (%) = 0.289037 (1.91719)

Particles: 124965 ave 126007 max 124347 min
Histogram: 4 0 0 0 2 0 0 0 1 1
Cells:      12500 ave 12500 max 12500 min
Histogram: 8 0 0 0 0 0 0 0 0 0
GhostCell: 1696 ave 1696 max 1696 min
Histogram: 8 0 0 0 0 0 0 0 0 0
EmptyCell: 0 ave 0 max 0 min
Histogram: 8 0 0 0 0 0 0 0 0 0
collide_modify      vremax 100 no
run 		    1000
Memory usage per proc in Mbytes:
  particles (ave,min,max) = 13.3896 13.3896 13.3896
  grid      (ave,min,max) = 3.73848 3.73848 3.73848
  surf      (ave,min,max) = 0.0601273 0.0601273 0.0601273
  total     (ave,min,max) = 17.2836 17.2836 17.2836
Step CPU Np Natt Ncoll Nscoll max 
    1000            0   999719      478      428        0           69 
    1100    1.5188251  1000362      482      434      119           78 
    1200    3.1278632  1000466      482      422      138           75 
    1300     4.769789  1000258      496      429      102           68 
    1400    6.2880991  1000850      486      439      119           74 
    1500    7.8050642  1001208      506      447      103           71 
    1600    9.4101341  1000756      477      411       92           67 
    1700    10.938024  1000834      470      428      116           68 
    1800    12.504573  1000050      503      443      103           70 
    1900    14.027482   999746      479      435      115           70 
    2000    15.597631   999970      502      451      128           80 
Loop time of 15.5977 on 8 procs for 1000 steps with 999970 particles

Particle moves    = 1003116823 (1B)
Cells touched     = 1120809552 (1.12B)
Particle comms    = 2907951 (2.91M)
Boundary collides = 0 (0K)
Boundary exits    = 2779034 (2.78M)
SurfColl checks   = 68389792 (68.4M)
SurfColl occurs   = 117996 (0.118M)
Collide attempts  = 15417194 (15.4M)
Collide occurs    = 11826282 (11.8M)
Particles stuck   = 0

Particle-moves/CPUsec/proc: 8.039e+06
Particle-moves/step: 1.00312e+06
Cell-touches/particle/step: 1.11733
Particle comm iterations/step: 1
Particle fraction communicated: 0.00289892
Particle fraction colliding with boundary: 0
Particle fraction exiting boundary: 0.0027704
Surface-checks/particle/step: 0.0681773
Surface-collisions/particle/step: 0.000117629
Collision-attempts/particle/step: 0.0153693
Collisions/particle/step: 0.0117895

Move  time (%) = 8.22284 (52.7184)
Coll  time (%) = 2.09558 (13.4352)
Sort  time (%) = 4.74273 (30.4067)
Comm  time (%) = 0.204996 (1.31427)
Outpt time (%) = 0.0418402 (0.268246)
Other time (%) = 0.289681 (1.85721)

Particles: 124996 ave 126052 max 124153 min
Histogram: 1 0 3 0 0 1 1 1 0 1
Cells:      12500 ave 12500 max 12500 min
Histogram: 8 0 0 0 0 0 0 0 0 0
GhostCell: 1696 ave 1696 max 1696 min
Histogram: 8 0 0 0 0 0 0 0 0 0
EmptyCell: 0 ave 0 max 0 min
Histogram: 8 0 0 0 0 0 0 0 0 0
